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Abstract 

We study in this paper the properties of a gas of fermions coupling to a 
U(l) gauge field at wavevectors q < A << kp at dimensions larger than one, 
where A << kp is a high momentum cutoff and kp is the fermi wave vector. 
In particular, we shall consider the e 2 — > oo limit where charge and current 
fluctuations at wave vectors q < A are forbidden. Within a bosonization ap- 
proximation, effective actions describing the low energy physics of the system 
are constructed, where we show that the system can be described as a fermion 
liquid formed by chargeless quasi-particles which has vanishing wavefunction 

overlap with the bare fermions in the system. 
PACS Numbers: 71.27.+a, 74.25.-q, 11.15.-q 
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I. INTRODUCTION 



In the last few years, there has been enormous interests in the study of U(l) gauge 
theories of fermionic systems in dimensions higher than one, as a result of appearance of 
effective gauge theories in the t — J type models [|l[] and also in the studies of the v = \ 
Fractional Quantum Hall state In particular, the behaviour of the systems in the 

strong-coupling regime where charge and current fluctuations are forbidden (confined) are 
of special interests. In a previous paper we studied the effects of a longitudinal gauge field 
(Coulomb interaction) on a system of spinless fermions in the limit e 2 — > oo. Within a 
bosonization approximation, we found that the ground state and low energy excitations of 
the system are described by a fermi liquid with chargeless quasi-particles [§J. In this paper 
we shall generalize our approach [|J to consider the e 2 — > oo limit of a system of spin S = 1/2, 
charge e fermions minimally coupled to a U(l) gauge field (A , A) in dimensions higher than 
one. The gauge field dynamics is described by an effective long distance action L gauge = F^ u 
with high momentum cutoff A << kp, where kp is the fermi wavevector. In this limit any 
non-uniform charge and current density fluctuations < p(q) >^ and < j(q) >^ with 
wave vector q < A in the system cost infinite energy and are forbidden. As a result any 
physical state \ip > that survives in this limit satisfies the constraints p{q)\ip >= and 
j{q)\ip >= 0, where p(q) and j(q) are the charge and current density operators, respectively. 

In this paper we shall show that by generalizing our previous approach 0], the projec- 
tion to the physical states \if) > where density and current fluctuations are forbidden can 
be carried out in the e 2 — > oo limit in a bosonization approximation |5|-f|, and physical 
(chargeless) single-particle operators that commute with both charge and current operators 
of the system can be constructed. The method we use is similar in spirit to the field theory 
scheme proposed by Shanjar and Murthy for the Fractional Quantum Hall effect, though 
the details are different. Within the approximation effective low energy Lagrangians that 
describe the chargeless particle-hole as well as single-particle excitations in the system are 
obtained. The physical sector of the system is described by a liquid of chargeless fermions 
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and corresponds to a kind of 'marginal' fermi liquid in the original fermion description. 
The organization of our paper is as follows: in section II we shall outline our mathematical 
formulation of the bosonization procedure where an approximate bosonized action of the 
system describing both spin and charge particle-hole excitations will be derived. In section 
III we shall study in detail the eigenstates and eigen-energies of the particle-hole excitation 
spectrum in the bosonization approximation for arbitrary values of coupling constant e 2 . In 
particular we shall show how the projection to the physical Hilbert space where charge and 
current fluctuations are forbidden is obtained in the e 2 — > oo limit. In section IV we shall 
consider the single particle excitations where S — 1/2, chargeless single fermion operators 
that commute with density and current operators are constructed and the equation of mo- 
tion for these fermion operators is derived. An effective action in terms of these chargeless 
fermion operators is constructed where we shall show that the system forms a marginal fermi 
liquid. Our results will be summarized in section V where some further comments will be 
given. 



II. MATHEMATICAL FORMULATION 

— * 

We shall work in the Coulomb gauge V.A = where the A component of the gauge 
field is integrated out, resulting in an instantaneous Coulomb interaction v(q) = 4ne 2 /q 2 at 
q < A between fermions. The resulting Hamiltonian of the system is, 

l r e - I 2 1 

H = S Y~ + 2I? ^ V (9)P($P(-<?) + ^photon, (1) 

where p ia and f ia are the momentum and position operators of the i th fermion with spin a. 
P(0) = ft +m Jk-q/2a is the density operator for the fermions. /(/ + )^'s are fermion 
annihilation (creation) operators. L d is the volume of the system and c is the velocity of 
light. 



Hphoton 

<f=^0,<?<A,A 



£ tv q (at(q)a x (q) + l), (2) 
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where a(a + )^\s are photon annihilation (creation) operators with wavevector q and polar- 
ization A. uj q = cq. We have set ft = 1 to simplify notation. Notice that the Hartree (q = 0) 
interaction energy does not appear in the Hamiltonian as in usual Coulomb gas problem 
where overall charge neutrality is maintained by presence of a uniform background of op- 
posite sign charges. After Fourier transforming and in second quantized form, we obtain 
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k 
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£ 2^ - c A ^ "> " m^M-A® A) (3a) 

la £ CT qX 

+ E(^)^A).l(-g,A), 
where j p (q) = J2k a (k / m ) f^ + ^ 2 a^k-q/2a * s ^ ne paramagnetic current operator, and 
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A(q, A) = f — I A)(a A (g) + a+(-g)), (3b) 



where £(<f, A) is the polarization vector of photons with momentum q and polarization mode 
A. Note that all sums over <f s are restricted to q ^ 0, |g| < A in Eq. (|3|) and in all 
the following equations. We have also replaced the fermion density operator p(r) by its 
expectation value Uq in the diamagnetic term in Eq. (|3"a|). This approximation can be 
justified in the e 2 — > oo limit where density fluctuations of fermions are forbidden. We shall 
discuss this more carefully in section IV. 

The difference between our model and usual Coulomb gas problem has to be emphasized 
here. In usual Coulomb gas problem the high momentum cutoff A is taken to be infinity, 
or satisfies A >> hp. In this limit a Wigner crystal where electrons position are fixed is 
expected to be formed (in 3D) when e 2 becomes large because the potential energy term 
becomes dominating at length scale ~ kp 1 . In our problem where A << kp, the interaction 
is effective only at length scale >> inter-particle spacing and formation of Wigner crystal is 
not warranted. In particular we found in Ref. [ [[|] that the ground state of the system has 
large degeneracies in the absence of the kinetic energy term and the system may remain in a 
liquid state when kinetic energy term is switched on. We note that bosonization method is a 



natural tool to study the liquid state of the problem in this limit [f7U8l,|10| . The bosonization 
procedure can be formulated most easily by introducing the Wigner function operators 
01 Pfc CT (^*) = ft+q/2afk-q/2a- 111 the path-integral formulation the Wigner operators are 
introduced in the Imaginary-time action of the system through Langrange multiplier fields, 
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dr 



f) k 2 

ha 



dr 2m 



ka,q 



+ JlT, v(q)p,:(q,T)pp(-q,T) - ^ ]T( LA ^ X ' T ' r) 



q,k,k' 



L d/2 



where /i is the chemical potential, 



a/ 



k,q\ 



2 

^ g 77,q 



r,/ 

^photon ' 



photon 



where S p hoton is the pure photon action, pj:(q) = -j^ H a P% a {(i) are density Wigner function 
operators, and A^ o .(g*)'s are Lagranger multiplier fields introduced to enforce the constraint 
that the Wigner operators are given by p% a (q) = f% +q -/2 a fk-cf/2cT- m particular, the original 
Hamiltonian ([!]) is recovered once the Xj: a {q) field is integrated out. For later convenience 
we shall also introduce spin Wigner operators cr^(q) = ■^J2a< J Pk a {q)- 
The photon action S' photon can be diagonalized easily to obtain 

SL^„ = V! i-ioJr,. + £l n )bt(q, ioL>„)b\(q, iojr,), (5a) 



'photon = (~ iuJ n + toq)K(& iUn)b X (q, IW„), 

q\,iu> n 



where 6(6 + )'s are new photon eigenmodes with eigenfrequencies VL q = JoJp + uj^, where 
ujp = \J Ann^e 2 / m is the fermion plasma frequency. The vector field A(q, A) can be written 
in terms of these new eigenmodes as, 



A(q,X) 



2tt 



m^)(h(q)+bt(-q)). 



(5b) 



Integrating out the fermion fields /(/ + )'s we obtain an action in terms of p% a {q) and 
X^(q) fields, 



S 







Trln 



GnX 



V2e ^ f k.A(q,\,iuj n ) ^ . 

Jl/2 2^ \ ~ )PkVli ~ lUJ n) + O photon , 



q,iu) n ,k,k' 



k,q\,iu) n 



where Fq is the free energy for an non-interacting Fermi gas. Go and A are infinite matrices 
in wave vector and frequency space, with matrix elements given by 
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G = $a<T>$k,k>go(k) g (k) 

rCCT, K (J 



(6b) 



and 



A 



, ,, , = S aa/ —=Xj: + p (k — k', iuj n — iuj n i) = 5 aa i—=\ k+k i (k — k 1 ), (6c) 



where k = (k,iu n ) and ^ = k 2 /2m — p. The Trln [1 — G \] term can be expanded in a 
power series of iX^ a (q) field, 



Trln 



1 - G a X 



—Tr 



GnX 



-Tr 



GqX 



2 1^ 
Tr 
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GnX 



+ 0(X 4 



Keeping terms to second order in A (Gaussian approximation), we obtain 



Trln 



l-G n X 



1 



2^ 9o{k + -)g (k - |)A feCT (g)A feCT (-g). 



(7) 



Notice that the first order term in A gives the usual Hartree self-energy and is excluded 
in our Hamiltonian. The iX^ a (q) fields in Action ( |6af ) can be integrated out in Gaussian 
approximation, resulting in an quadratic action S in terms of p%(q), and photon fields 

only. We obtain S = S p ^ p + S a , where 



k,k',q,iuj 

\[2e 



2v(q) 
L d 



p%(q, iu n )pp(-q, -iu) n ) 



(8a) 



,k.A(q,X,iuj n ) . , 

L d/2 2- ( )Pk(Q>- iaj n)+Sp h oton 



k,q\,iij n 



m 



and 



S * = -l 7^-- — ^k(^ iu n)^(-q,-iuj n ), 

k,q,tu} n u/t 



(8b) 



where 



1 



fc— g/2 n k+q/2 



X k(q, = j- 9o{k + q/2, iuj n + iVL n )g Q (k - q/2, ifl n ) 

P in n iuj n - -fr 



(8c) 



rif: = 0(—^) at zero temperature is the free fermion occupation number. S p ^ p and S a can 
be expressed in terms of canonical boson fields by introducing 

Pk(Q,^n) = yj\&n(q)\ (e(A^)b+(q,icu n ) + e(-A^)b R (-q, -iu n )) , (9a) 
a % (q, iu n ) = y/\A^(q)\ (e{A^(q))si{q, iuj n ) + 6{-A u (q)) Sj ;{-q, -iu r , 

where A^(q) = n%_-, 2 — n^ + ^, T Correspondingly, we also have 



Pi (-q, -iuj n ) = J\A^q)\ (6{A % {$)b % {q,iuj n ) + 9(-A^)bt(-q, -iu n )) , (9b) 



a^-q, -iu n ) = ^JlA^qJl ^(A^(g))s^(g, iu n ) + e(-Aj;(q))s^(-q, -iu n )) . 

Putting eqs. @ back into S p - P , we obtain after some straightforward manipulations, 

1 Ifc.gl 1 
S p-p=2 1 (- iuJ n H — )b%{q,iu n )b%{q,iu n ) + - (-i^n + Q q )bt(q,iuj n )b\(q,iuj n ) (10a) 

k,q,iut n &,i<*>n 

+J1 II v(q)^\A^A^\9(A^)e(A p m x (b+(q,iu n )b p (q,iu; n ) 

k,k',q,iu> n 

+bi(q,iu} n )b+-,(-& ~ iuJ n) + b_i;(-q, -iu n )bp(q,iuj n ) + b_^(-q, -iuj n )b + ^,{-^ -i^n)) 

k,q\,iuJ n 

x(b x (q,iuj n ) +b^(-q,-iu n )), 



and 



1 \k.q\ 
S ° = o (-iuJ n + -^)s^(q,iuJn)s l :(q J iuJn), (10b) 



k,q,iu) n 



where 



m^qA) = —(^-)^k.i{q,X). 
m \l q 

S p ^ p describes a coupled system of (fermion) density particle-hole pair excitations and pho- 
tons, whereas 5*0- describes a system of free (fermion) spin particle-hole excitations. The 
density- and spin- particle-hole pair excitations are described by boson fields b{b + )^(q) 
and s(s + )^(q), respectively, satisfying usual boson commutation rules [b(s)j^(q), b(s)i(q')] = 



SktfScffl and [b(s) a ,b(s)p] = [b(s)^,b(s)t] = 0, etc. In this form the dynamics of the original 
fermion system is described completely in terms of charge- and spin- boson fields(bosonized). 

Notice that we have restricted ourselves to the Gaussian approximation in deriving S p - P 
and 5*0-. Higher order interaction terms between bosons will appear in a cumulant expansion 
of the Ar(g) fields [|It]]. The convergence of the cumulant expansion is formally controlled 
in our model by the small parameter e ~ (A/ hp), where a term of order [6 m 6+( ?1_m )] in the 
cumulant expansion has to leading order a term ~ e 2 (d-i)x(«-i) j n eac h increasing order of 
cumulant expansion |lT|fl. Notice however that the smallness of e does not guarantee that 
all calculated physical quantities will converge uniformly in the cumulant expansion. In 
particular, we shall see that infra-red divergence appears in higher-order expansion of the 



transverse gauge field [12|,[T3j . We shall discuss in more detail the consequences of the infra- 
red divergences in section V. Notice also that in the g* — ^ limit, A^(g) — > — S(e^— //)(—) and 
the usual "tomographic" bosonization procedure based on subdivision of Fermi surface into 
disjoint patches at small q is recovered [§-§]. Our bosonization procedure can be viewed as a 
generalization of the tomographic bosonization method for small wave vector q to arbitrary 
values of \q\ < A. 

It is straightforward to show that the Gaussian approximation is essentially the same 
as usual RPA theory for interacting fermions @J6|-[8| , except the additional assumption that 
particle- hole excitations can be treated as independent bosons in bosonization theory [|J. 
Note that particle-hole pairs are not all independent in a fermion system because of the Pauli 
exclusion principle. The assumption of independent bosons is a major approximation in the 
Gaussian theory in two- and three- dimensions. This is in contrast to what happens in one 
dimension where the entire particle-hole excitation spectrum can be represented rigorously 
by bosons || when the fermion spectrum is linearized near the Fermi surface. 

Despite the approximations in the Gaussian theory, the bosonized form of the action 
has the advantage that within the approximation the full excitation spectrum as well as the 
wavefunctions of the system can be obtained easily. This allows us to study the properties 
of the system in great detail, as we shall see in the following. 
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III. DENSITY PARTICLE-HOLE EXCITATION SPECTRUM 



The eigenstates and eigenvalue spectrum described by the action S p ^ p can be obtained 
by diagonalizing the bosonized action ( |10a|) using a generalized Bogoliubov transformation. 



We introduce for each wave vector q the Bogoliubov transformation |14 



h(q) = Y [ a kk'lk'{q) + Pkk'J-k>(-q) 
b-k(-q) = Y [ a kk>l-k'(-q) + Pkk'lk'io) 



(11a) 



where k = k,X(k' = k',X'), and with correspondingly, 

lk(q) = Y [ a k*kh'(q) ~ Pk>k b -k'(-q) 
fe' 

i-k{-q) = Y [ a k>l b ~k>(-q) - PvkK'io) 

where we require that the 7(7 + )fc(<?) operators diagonize the Hamiltonian, i.e. 

H p _ p = Y E Mlt{$lM + E G , 



(lib) 



(12) 



where Ek(q) are the eigen-energies and Eg is the ground-state energy of the system. The 
eigenstates may represent dressed particle-hole excitations (k = k) or dressed photons (k = 
A). Notice that additional collective modes may appear in the system and is also included 
in the sum J2k- The matrix elements a and j3 satisfies the orthonormality condition 



>«) >«)* _ o>«)o>(<)« 

a kk" a k'k" Pkk" Pk'k" 



Yj l a kk"Pk'k" Pkk va k'k" 



Skk', 
0. 



(13) 



Writing down the equation of motions for bkiffls in terms of 7(7 + )fc(<?)'s we obtain the 
Bogoliubov equations 



(E k ,(q) )a T 



m 



kk' 



>(<) , Q <>>^ 

\k' ) 



(14a) 



I Ml. 



(E k ,(q) + J ^)/3gi >) = -0(Ar(gW|Ae(g)| 



m ' kk 



^ 37 |l^ mE (,-,A)(aJ'i, <) - (14b) 



4 L d 

k" 



and 

(^(?) - njaJi? = : ^E^(9'^)«( A if-(9))>/l A iP'(9)K a ^ ) -4^)' ( 14c ) 

fc" 

(Mfl + ^)^ >} = ^ A)6>(A*(fl)^|A*(fl|(a><<>' - /?<£>). (14d) 

fc" 

In general we find that there exists two kinds of solutions to these equations: (i)particle- 

— # 

hole continuum, with E^(q) = \k.q\/m, and (ii) collective modes, including (renormal- 
ized) photons with energy E\(q) determined by the eigenvalue equation E\(q) 2 — Q 2 — 
^rXt(q, E\(q)) = 0, and plasmons, with energy E (q) given by the eigenvalue equation 
1 - v(q)xo(q, E (q)) = 0, where 

2 



= £( fc *-**)x fc(g>)» 



is the paramagnetic transverse current susceptibility, k t = k — q(k.q) and 

2 

Xo(5>) = ^dZ)Xofc(9' w )- 

k 

is the Lindhard function. The energies of these collective modes are outside the particle-hole 
continuum. We now consider the solutions in more detail. First we consider the particle-hole 
excitations. After some lengthy algebras, we obtain, 



g(Ag(g))g(A p (g))^/|Ag(g)Ag y (g)| 



V m m / 



^Xk' 



(Ik'.gVm + ^)t(^)^(A g ,(g))^A^ . 

L^((|fc'.gi/^) 2 - ^ - \k'4/m) [Q) ' th 
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for the particle-hole continuum spectrum k', where v e tAq,uj) = v(q)/(l — v(q)xo(q,u) is the 



RPA effective interaction, and A eff (q,k') = (^k' t ) /{{\k' .q\/m) 2 - Q 2 q - ^f Xt (q,\k' .q\/m)) 
We also obtain ar(/3)< & = -a(/%„ and = a(/3)s> 



For the collective modes, we obtain for the longitudinal (plasmon) mode, 

^ i 2e(A~M)J\A~M\ 



a 



/,■(") 



<* 



L d/2 



{EM - 



[Ml - 



9oj 



(15b) 



L d/2 



«A0 = A 



<* 
AO 



(E (q) + 
0. 



1Mb 



dxo0?>^) 



<9w 



and «(/3)| Q = a(/3)^ Q , a(/3)^ = 0. Notice that in bosonization theory, plasmon arises from 
non-perturbative interaction effect and cannot be obtained from analytic continuation of 
the non-interacting boson modes. Similarly, we also obtain for the (renormalized) photon 
modes, 



a 



kX 



<* 
kX 



e(A^q))^\^M\ kUq)C(q) \ 

EM + ^ \ Ld/2 ) 

i 

5\A' < 



7T 



(E\(q) - Vt q )m \n q 
e I tt 



Xt(q,E x (q))C(q), 



S 



XX' 



ftp f-*,n\ \7r) Xt(q,E x (q))C(q), 



where 



C{q) 



' 4?re 2 \ r dxt{q>u) l 



(15c) 



and with a(f3)f, x = a(J3)% x and «(/3)| A = -a(fi)g x - 

Next we examine the solutions of the bosonized Hamiltonian in the e 2 — > oo limit. First 

2 

we consider the collective modes. Using the result that Xo(q, — > an d Xt(q,u) — > 



2n q 2 
dw 2 



< e > in the limit uj » kpq/m, where hq is the fermion density [pjj and < e > is the 
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average kinetic energy per fermions in the free fermion ground state, it is easy to see that 
in the limit e 2 — > oo, the collective mode frequencies E (q), E\(q) all approaches the plasma 
frequency uip = ( 47r "° 6 )i. Notice that up — ► oo as e 2 — > oo, indicating that the collective 
modes are outside the physical spectrum in this limit. 

Despite the vanishing of collective excitations in the physical spectrum, the particle- 
hole excitation spectrum with excitation energies \k.q\jm survives in bosonization the- 
ory in the limit e 2 — > oo. In this limit v e ff(q,\k'.q\/m) — > — l/xo(<7, |&'-<7|/ m ) an d 
k.A e ff(q,k') — > — 2k t .k' t / {mriQ + Xt(9> l^'-9l/ m ))- The coefficients agg/'s and /%£/'s remain 
regular, indicating that the particle-hole excitation spectrum survives under effect of con- 
finement. It is instructive to examine the charge and (transverse) current fluctuations car- 
ried by the (eigen)-particle-hole excitations by examining the commutators \p(q), 7g(<f )] and 
[Jt(<7%7fc(<f)], wh ere 

1M = E —piM - (L d/2 ) E —A-Z a), 

- m ^ m 

fccr A 

is the transverse current operator. In particular, we expect that these commutators should 
vanish in the e 2 — > oo limit, where charge and current fluctuations are forbidden. Using Eqs. 
(0), ( [Lib ), ( |15a| ) and usual boson commutation rules, it is straightforward to show that 



— e 



and 



1 - v{q)Xo{q, \k.q\/m) 



&®nm = Sv^8(±M\/\±Mx ( lT* 1? j^^r I • (17) 



m 

both vanishes in the limit e 2 — > oo. 

Let us summarize our results obtained so far from bosonization theory. Within the Gaus- 
sian approximation, we obtain a RPA-like charge excitation spectrum with both collective 
modes and particle-hole excitations, and a free-fermion type spin particle-hole excitation 
spectrum. As the coupling constant e 2 increases, the energies of the charge collective modes 
rise continuously to infinity whereas both the charge and spin particle-hole excitation spec- 
trum remains intact. The charge and current fluctuations carried by the charge particle- hole 

12 



excitations are projected out gradually as e 2 increases, resulting in chargeless particle-hole 
excitations in the confinement limit e 2 — > oo. Notice that within the Gaussian approxima- 
tion, the confinement state analytically continues to the usual Fermi liquid state and there is 
no phase transition in between. The theory thus suggests that the confinement state of a gas 
of fermions is a fermion liquid state formed by spin S = 1/2 and c/iarge/ess-quasi-particles. 
It also suggests that this is a rather unusual fermion liquid state, since bare fermions in the 
system carries spin 1/2 and charge e, and the quasi-particles must have vanishing overlap 
with bare fermions if they carry zero charge. The nature of the liquid of chargeless fermions 
will be examined in more details in next section, where we shall construct explicitly "charge- 
less" fermion operators that describe the dynamics of the system in the limit e 2 — > oo. The 
limitations of Gaussian theory will be uncovered in the process. 



IV. SINGLE-PARTICLE PROPERTIES AND LOW-ENERGY EFFECTIVE 

LAGRANGIAN 

In bosonization theory for one-dimensional systems, the single particle properties of the 
system can be determined once a rigorous representation of the single-particle operator in 
terms of density operators Pl(<?) an d Pn(q) are obtained ||. In higher dimensions, this 
procedure becomes inadequate for two reasons: (l)The corresponding procedure requires 
that the bosonized representation of single-particle operator ipbo{r) satisfies the commutation 
relations 



&b*(r),Pk(q)\ = e-^ 12 ^ J d d r' e^' 2 ^ 



W(0 



for all possible momenta k and q. We have not been able to find a representation which 
satisfies this criteria 0,§,|II],|J , and even if we can find such an representation, the theory 
would still be approximate because in dimensions higher than one, the boson representation 
using Wigner operators is not exact, (2)more importantly, unlike in one dimension where 
the elementary excitations are collective density- and spin- wave modes, bosonization theory 
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suggests that in dimensions higher than one the particle-hole excitation spectrum is fermi- 
liquid like, implying that fermionic quasi-particles exist in dimensions higher than one. It is 
important to construct directly quasi-particle operators in this case. 

To construct the quasi-particle operators we start from the equation of motion of the 
bare fermion operator ip a {f) = J2k e%h ' r fk a a ^ imaginary time, 



Or 



[#,W(rO] (18) 
1 1 - - ie 

q q\ 



where we have replaced the fermion-density operator p(r) by its expectation value n in 
the diamagnetic term as we have done in sections II and III in deriving ([Tj|). We shall 
discuss the validity of this approximation at the end of this section. In the bosonization 
approximation, the operators p(q) and A(q, A) can be decomposed as p(q) = p p h(q) + Pc(q), 
and A(q, A) = A p h{q, A) + A c (q, A), where v(q)p p h(q) and eAph ^' x ^ describes the coupling 
of the fermion to particle-hole excitations through the scalar and vector field fluctuations, 
respectively, whereas v(q)p c (q) and eAc ^^ describes coupling of the fermion to collective 
modes (plasmons and photons). It is straightforward to obtain 



>(q) Pph (q) = j:j\ A MmM 



m 



m 



(19a) 



,(q)p c (q) = 



du 



u)=E (q) 



7o"(5) +7o(-5) 



and 

eA ph (q, A) 



in 



V2L d / 2 



In particular, it is easy to see that in the e 2 — > oo limit, the interaction between fermion 
and particle-hole excitations in both longitudinal (v(q)p p h(q)) and transverse (eA p h(q, A)/c) 
channels remains regular and finite, and divergences in the e 2 — > oo limit appear only through 
the interactions between fermion and collective excitations. Notice that infra-red divergences 
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in the interaction between fermion and particle- hole excitations also exist. However, they 
are not results of taking e 2 —>■ oo and are not considered here. We shall discuss the effects 
of infra-red divergences in the last section. 

The divergence in interaction between fermions and collective modes suggests that we 
have to eliminate these interactions first to study dynamics of (physical) fermion operators 
in the limit e 2 — > oo. To eliminate these interactions we look for a canonical transformation 
for the fermion operator [|T5 1 , 



Vw(r) =e* f tya(r)- (20) 

with corresponding equation of motion, 

[if, $ Pa {f)\ = [H, e^]^(r) + e^[H, Vv(f)], (21) 

where we shall choose 0(f) such that [H, e^} cancels the interaction term between fermion 
and collective modes, i.e. 



[H, e^]Vv(r) =e* w 



Td E v(q) P Me- iq - r M?) + ^Jf, E e n - r MQ, A)-V^(r) 

q qX 



Furthermore, we shall assume that 0(f) depends linearly on the collective mode operators 
7o(A)(<?) and 7 ( ^ A - ) (g)'s, and that H can be replaced by the bosonized Hamiltonian ( |T2"D in 
evaluating the commutator [H, e^]. Notice that the last assumption is justified in the e 2 — > oo 
or q — > limit, where both the longitudinal and transverse collective mode excitations 
become exact. 

With these conditions, it is straightfoward to show that 0(f) = 4>i(r) + W(f).V, where 

i p-iq-f 

m = —2 E r - ; i (7o + (fl - 7o(-9)), (22) 



L d l 2 . 

i E (q) 



where 0; and describes the dressing of fermion by (longitudinal) plasmon modes and 
(transverse) photon modes, respectively. Notice that [<f>i(r), W(r > ).V] = 0, because of de- 
coupling between transverse and longitudinal fluctuations. 

15 



To show that ipp a {r) describes chargeless fermions in the e 2 — * oo limit, we examine the 
commutation relations between ipp a {f) and density and (transverse) current operators. It is 
straightforward to show that 

/ 



[ep(q),ip P<T (f)] = e 

and 



Eo(q) 



2Xo(g ' Eoiql) l|e**lM*0, (23a) 



du> 



uj=E (q) 



iJt(q), lMr)l = ^ (l " IxM E x {q)) + «om]) e^VrfMO- (23b) 

where Vt = V — <?(<?. V). It is easy to see that both commutators vanish in the limit 
e 2 — > oo, when _E (a)(<?) — ► — ► oo. It is also easy to see that the ipp a operator has 
the same commutation relation with spin density operator a(q) = JZu^tio) as bare fermion 
operator ip a . These results together imply that the dressed single particle operators ippaif) 1 ^ 
represent spin S — 1/2, 'chargeless' fermions in the limit e 2 — > oo. 

To show that ipp a {r) and xjj Pa [r > ) represent independent physical fermionic excitations in 
the system when f ^ r* we examine the commutation relation between the dressed fermion 
operators. We obtain 

in the limits e 2 — > oo and \ f — f*| — > oo, where 

6 = (^4). te(^)e^(V^(r-))e^ + ^) - (V^(0)e^(rV +( ^ . 
|r — r | L J 

The main contribution to the commutator comes from nonzero commutation relation be- 
tween the plasmon cloud around one chargeless fermion (4>i(r)) and the other bare fermion 
(■00-/ (f*)). The non- vanishing commutator between chargeless fermions separated by large 
distance reflects the nonlocal nature of the chargeless fermion operator. Fortunately, the 
power law decay of the commutator between ipp a operators separated by large distances at 
dimensions larger than one indicates that they can be used as starting point to construct 
independent single fermion operators when describing the dynamics of the system at long 
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distances. Notice that in one dimension, such a construction is not possible because of dis- 
tance independence of the commutation relation. In fact, the only fermionic operators that 
commute with density operator p(q)'s are the ladder operators || which raise or lower the 
number of particles in the system by one. There are only four independent ladder operators 
in the system, which change the number of left-going and right-going spin-cr fermions and 
cannot be used to construct local quasi-particle excitations. 

The relation between the fermion liquid formed by the chargeless and bare fermions 
can be seen by examining the relation between the two ground state expectation values 
< ip^(r)tl> a (r') > and < ipp a {r)ipp a {f) >. Assuming that the collective modes can be 
treated as independent excitations, we obtain to leading order in the projected Hilbert 
space, 

< ipj (r)ipo.(r) >~ e n ° ' Jl - 1 ™ wv < ipp a (r)ipp a (r) > . (24a) 

The singular behaviour in the exponential factor comes from plasmon contributions and 
photons do not contribute to this order. Notice that the average density of bare fermions 
and dressed fermions are the same and the fermi surface volume of the dressed fermions is 
exactly the same as that of the bare fermions. It is also easy to see that 

< tpt^ip^r) >~ e " (S ^" )3( -^- 1 ddq ~? ] < ip+(r)ip a (r) >, (24b) 

which vanishes in the e — > oo limit, indicating that the bare and chargeless fermions have 
zero wavefunction overlap as expected. Using Eq. ( |24a| ) , it is straightforward to show that 
the bare-fermion occupation number n(k) has no discontinuity across fermi surface in the 
e 2 — > oo limit (3D), even if the chargeless fermions form a fermi liquid, in agreement with 
the marginal fermi liquid picture. 

The bosonization result we obtained in previous section suggests that in the limit e 2 — > oo 
the low energy physics of our system may be described as a fermi liquid of the chargeless 
fermions. To see whether this is the case we consider the equation of motion for the chargeless 
fermion operators ipp a , assuming that they can be treated as canonical fermions. It is 
straightforward to show that 
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= i ( V2 ^W " [V 2 0(r)- : (V0(f)) 2 #P CT (r) - 2V0(r).V^P CT (r)) (25) 
-^E«(ff)ft*(ff)e**Vp ff (T=) - ^£^"^(9, A).V^(tO, 

where we have neglected constant energy terms coming from normal ordering of operators 
inderiving Eq. (|25|). Notice that direct couplings between fermions and collective modes 
are absent in the equation of motion of ipp a (r). However, interactions between chargeless 
fermions and particle-hole excitations remain in the equation of motion. Moreover, an 
indirect coupling to collective modes is also generated from the fermion kinetic energy term 
as is in the similar "small polaron" problem \IE\. It is easy to see by direct power counting 



of e 2 in the cj){r) field that the coupling of the "dressed" fermion to collective excitations 
through the kinetic energy term is much weaker than the original fermion-plasmon and 
fermion-photon couplings. In particular, the self-energy correction of chargeless fermions 
from </>(r) fields remains finite in the e 2 — > oo limit. 

We shall now study the dynamics of the chargeless fermions in more detail. First we 
write down an effective action for the chargeless fermion operator ipp a {r), the effective action 



is constructed with the requirement that it reproduces the equation of motion fl25|) for i/jp a . 
It is straightforward to show that the correct action is 

d V 2 



rP 



dr 2m V + iA eff (x,T).V + (j) eff (x,T) 



ip Pa (x, r), 
(26a) 



where (j) e ff(q) ~ v (<l)Pph{.<i) and A e ff(q) ~ eA p h(q)/m, with dynamics given by 

S e ff((f>) = -i^Xo(g,^)|0 e //(<?,^)| 2 , (26b) 
S eff{A) = --J2Xt(q,iu)A ef f(q,u).A eff (-q,-iLj) 

where Xt{q, iu) — Xt{q, ioo)+mno is the total transverse current susceptibility of free fermions. 
We have neglected in S e ff(ipp, ipp) the collective mode contributions to the dynamics of ipp a 
through the kinetic energy term. The residue plasmon-fermion interaction induces a weak 
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short-ranged effective attractive interaction between chargeless fermions and may lead to 
superconductivity. We shall not discuss this possibility in this paper. 

Fermions interacting with gauge fields with effective interaction ( pE| ) have been studied 
in detail by a number of authors [3,|3i |12| , |13| and it is believed that interaction between 



fermions and transverse gauge field with effective dynamics ( |2B| ) may lead to break down 
of fermi liquid theory and the formation of a marginal fermi liquid. If this is the case, the 
ground state formed by the chargeless fermions would itself be in a marginal fermi liquid 
state and our system can be viewed as a "double" marginal fermi liquid state of the original 
fermions. This is in contrast to (Gaussian) bosonization theory which predicts that the 
chargeless fermions should form a fermi liquid state. The failure of Gaussian theory in 
producing a marginal fermi liquid state shows the limitation of Gaussian theory. This point 
will be elaborated further in the last section. 

It is interesting to make comparison between our approach and the usual treatment of 
the U(l) gauge field in t-J model where the gauge field effectively imposes the constraints 
Ps{q) + Ph{q) = and j s (q) +j h (q) = |, where p s (J s ) and p h (J h ) are the spin and hole 
densities (current densities) operators, respectively. In the usual treatment [0,0], the spins 
and holes were integrated out at Gaussian level, resulting in an effective action for gauge 
fields very similar to Eq. (|26| ) except that the longitudinal and transverse response functions 
Xo(g, uj) and Xt(q, are replaced respectively by the sum of the longitudinal and transverse 
response functions from the spin- and hole- components, respectively. It is then assumed 
that the spinons and holons in the projected Hilbert space can be treated as free particles 
interacting with the gauge fields with dynamics given by S e ff [p]]. Although our analysis 
here includes only fermions and is not directly applicable to the t-J model, the form of our 
low energy effective Lagrangian is in agreement with this assumption, suggesting that the 
spinons in the usual treatment of t-J model are strongly related to our chargeless fermions. 
A detailed treatment of the t-J model using our approach will be reported in a different 
paper. 

Finally we examine the validity of replacing the fermion density operator p(r) by the 
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average fermion density no in the diamagnetic term when we derive the equation of motions 
for the bare and chargeless fermions. We shall now show that this approximation is justified 
in the e 2 — > oo limit for the chargeless fermions tpp a - We look at the correction term in the 
equation of motion for chargeless fermions, [AH,i/)p a (f)], where 



AH = 




It is easy to show from Eqs. (|19a|) , (|23a|) and (|2~2"D that [p(q), ipp a {r)] ~ ■^i'Pa, eA/m ~ 
c 1//2 (7a + 7a) an d [ e ^4(^)> ^Paif)] ~ V^Po- in the e 2 — > oo limit. As a result we have 

[AH,^ Pu (r)] ~ {A{?).V1> Ptr {f)){p{f) -n ) + 0(~). 

The first term vanishes in the e 2 — > oo limit if the matrix elements < n\(p{f) — n )|m > 
vanish faster than 1/e 1 ' 2 in the physical Hilbert space of interests. This will be the case if the 
physical Hilbert space is spanned by the chargeless fermion operators we construct, because 
of the commutation rule between density operator and chargeless fermions (Eq. ( |23a| )). 

V. SUMMARY 

Using a bosonization approach we studied in this paper a gas of spin S = 1/2 fermions 
interacting with a U(l) gauge field with high momentum cutoff q < A << kp. In particular 
we consider the e 2 — > oo limit where the gauge field becomes confining. We have analysed 
the problem in two steps. Within the Gaussian approximation we show in section III that 
a liquid state solution is obtained where the particle-hole excitation spectrum of the system 
is always fermi-liquid like, with the charge carried by the particle-hole excitation vanishing 
continuously in the e 2 — > oo limit. To verify the above picture we construct in section IV S = 
1/2, chargeless fermionic operators that can be used as the starting point for constructing 
quasi-particles in the system. We find that the dynamics of these "physical" single-particle 
operators are governed by a Lagrangian very similar to the effective Lagrangain obtained 
in conventional treatment of gauge field model |TJ. The solution we obtained describes a 
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kind of 'double- marginal' fermi liquid state where (i)the chargeless fermions have vanishing 
wavefunction overlap with bare fermions in the system and (ii)they form a marginal fermi 
liquid state themselves. 

It is useful to distinguish the effects of longitudinal and transverse gauge fields on our sys- 
tem, in particular how they contribute in the e 2 — > oo limit, to the construction of chargeless 
fermions and their effective dynamics. First we consider the longitudinal gauge field. The 
main effect of longitudinal gauge field appears through fermion-plasmon interaction which 
is singular in the limit e 2 — > oo. In dimension d = 2, there exists an additional infra-red 
singularity associated with the response of plasmon field (orthogonality catastrophe effect) 
when charged particles are added to the system |17j . The effect exists for any finite values of 



e 2 > 0. The interaction between fermions and particle- hole excitations is finite and regular 
for any value of e 2 . The singular fermion-plasmon interaction and the infra-red singlarity at 
dimension d = 2 are eliminated together with our canonical transformation fl2~0|) ||], and 
the remaining interaction between the chargeless fermions and longitudinal gauge fields is 
regular and introduces no further singlarities in the chargless fermion system ||. 

The situation is however, quite different with the transverse gauge field. There exists no 
infra-red singularity in the fermion-photon interaction at both two- and three- dimensions 
when e 2 is finite, and the singularity associated with the fermion-photon interaction in the 
limit e 2 — ^ oo is in fact weaker than the corresponding singularity associated with fermion- 
plasmon interaction, as can be seen from direct power counting of e 2 in <pi and W fields. 
As a result, the canonical transformation ( p0|) is dominated by effect of plasmons in the 
limit e 2 — > oo. On the other hand, the interaction between bare fermions and particle- hole 
excitations through transverse gauge field is infra-red singular for any values of e 2 in both 
two- and three- dimenions and the singularity carries over to the effective interaction 
( PSD between dressed fermions and transverse gauge field. It is believed that the infra-red 
singularity in interaction between fermions and particle-hole excitations through transverse 
gauge field will lead to the formation of a marginal fermi liquid state 0J^,|T^,|T^ . 



The possibility of formation of an marginal fermi liquid state of the chargeless fermions 
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points out the limitation of bosonization theory in Gaussian approximation, where a fermi- 
liquid type particle- hole excitation spectrum is always obtained in the limit e 2 — > oo. It 
is clear that higher-order terms in cumulant expansion of transverse gauge field has to be 
included in bosonization theory to obtain the correct marginal fermi liquid behavior [0, 
and a complete theory where particle-hole and single-particle excitations are treated self- 
consistently is still missing. 

An immediate question that comes with this observation is that we have used the exci- 
tation spectrum in Gaussian approximation to construct chargeless fermion operators in our 
paper. If the particle-hole excitation spectrum is incorrect, our results based on chargeless 
fermion operators ipp a become questionable! Fortunately although the particle-hole excita- 
tion spectrum in Gaussian theory is questionable, the plasmon and photon excitations we 
obtained are exact in the q — > or e 2 — > oo limits. In particular our canonical transfor- 
mation p0| ) for the chargeless fermions ipp a involves only the collective mode operators 
and is independent of the Gaussian approximation! However, it has to be cautioned that in 
deriving the effective action S e ff(ipp,ipp) which describes scattering of chargeless fermions 
with particle-hole excitations, it is assumed that the eigen-particle-hole excitations are de- 
scribed correctly by Gaussian theory. Strong modifications may occur if the spectrum of 
particle-hole excitation is modified strongly in the correct theory. 

Another more fundamental question is how valid it to treat our chargeless fermions 
as canonical fermions, given that they obey a rather non-trival commutation relation. In 
particular, it is possible that the unusual commutation relation between fermions may lead 
to additional non-fermi liquid behaviour. At present we have no answer to these questions. 

This work is supported by HKUGC through RGC grant HKUST6124/98P. 
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